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Abstract 

We give a short review of the recent development in the investigation of the vacuum 
expectation values of the bulk and surface energy-momentum tensors generated by quantum 
fluctuations of a massive scalar field with general curvature coupling parameter subject to 
Robin boundary conditions on two codimension one parallel branes located on (D + 1)- 
dimensional anti-de Sitter (AdS) bulk. An application to the Randall-Sundrum braneworld 
with arbitrary mass terms on the branes is discussed. 


PACS numbers: 04.62.+V, ll.10.Kk, 04.50.+h 


1 Introduction 

Recent proposals of large extra dimensions use the concept of brane as a sub-manifold embed¬ 
ded in a higher dimensional spacetime, on which the Standard Model particles are confined. 
Braneworlds naturally appear in string/M-theory context and provide a novel setting for dis¬ 
cussing phenomenological and cosmological issues related to extra dimensions (for reviews in 
braneworld gravity and cosmology see Refs. [J). The investigation of quantum effects in 
braneworld models is of considerable phenomenological interest, both in particle physics and 
in cosmology. The braneworld corresponds to a manifold with boundaries and all fields which 
propagate in the bulk will give Casimir-type contributions to the vacuum energy, and as a result 
to the vacuum forces acting on the branes. In dependence of the type of a field and boundary 
conditions imposed, these forces can either stabilize or destabilize the braneworld. In addition, 
the Casimir energy gives a contribution to both the brane and bulk cosmological constants and, 
hence, has to be taken into account in the self-consistent formulation of the braneworld dy¬ 
namics. Motivated by these, the role of quantum effects in braneworld scenarios has received a 
great deal of attention (see, for instance, references given in (2j). In particular, the investigation 
of local physical characteristics, such as expectation values of the bilinear field products and 
the energy-momentum tensor, is of considerable interest. In addition to describing the physical 
structure of the quantum field at a given point, the energy-momentum tensor acts as the source 
in the Einstein equations and therefore plays an important role in modelling a self-consistent 
dynamics involving the gravitational field. In this paper based on Refs. HE], we give a short 
review of the recent progress in investigations of the vacuum expectation values (VEVs) for 
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both bulk and surface energy-momentum tensors in braneworlds with AdS bulk (for the vacuum 
energy-momentum tensor in the Randall-Sundrum braneworld see also 21). The corresponding 
results for the geometry of two codimension one parallel branes located on (D + l)-dimensional 
background spacetime AdSox+i x £ with a warped internal space £ are presented in [jjJJ. 


2 Bulk energy-momentum tensor and vacuum interaction forces 


Consider a massive scalar field (p(x) on background of a (D + 1 (-dimensional AdS spacetime 
(AdSp+i) with the line element 

ds 2 = gikdx l dx k = e~ 2kL>y r] tiI/ dx fi dx l/ — dy 2 , ( 1 ) 

and AdS radius given by 1/ko- Here r/ M , y = diag(l, — 1 , ..., — 1 ) is the metric for the D- 
dimensional Minkowski spacetime, i,k = 0,1,..., D, and /j, u = 0,1,..., D — 1. The corre¬ 
sponding field equation reads 

([ViV* + m 2 - (D(D + 1 )k 2 D \ <p(x) = 0, (2) 

where the symbol V,; is the operator for the covariant derivative, £ is the curvature coupling 
parameter. We will assume that the field obeys Robin boundary conditions on two parallel 
infinite plane boundaries (branes), located at y = a and y = b, a < b: 

Ay + Bydy^j <f(x) = 0, y = a, 6, (3) 

with constant coefficients A y , B y . The eigenfunctions of the problem are decomposed into 
standard plane waves with the wave vector k in directions parallel to the branes and into the 
transverse part which is a cylindrical function of the argument uz, where 2 = e k ° y /kr>. The 
eigenvalues u = u u ^ n , n = 1 , 2, ..., are solutions to the equation 

gl ab) ( uz a ,uz b ) = 4°) (uz a )Yj b) ( uz b ) - yj a) ( uz a )J {b) ( uz b ) = 0, (4) 

where J v (x), Y v (x) are the Bessel and Neumann functions of the order 


v = \J (D/2) 2 - D(D + 1)C + m 2 /k 2 D . (5) 

In Zj = e k °i / kn, j = a, b, and the notation 

F®(x) = AjF(x) + BjxF'ix ), ( 6 ) 

is introduced for a given function F(x), with the coefficients 

Aj = Aj + Bjk[)D/ 2, Bj = Bjkr>. (7) 


The VEVs of the physical quantities bilinear in the field operator (Wightman function, field 
square, energy-momentum tensor) are expressed in terms of the mode-sums which contain the 
summation over the eigenvalues u v ^ n . The application of the generalized Abel-Plana formula 
0 to this sum allows us to extract the boundary-free AdS part and to present the boundary 
induced parts in terms of exponentially convergent integrals. In particular, the VEV of the 
energy-momentum tensor in the region between the branes is presented in the form 


(0\T k \0) = (0 s \T k \0 s ) + (T k )^- 


k° +l z D 5 k 


2 D ~ 1 ir^r (f) Jo 


r 


duu D 




(' uz a ,uz b )F W G ( J\ 


UZ n ,UZ 


( 8 ) 
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where j = a,b provide two alternative representations. In this formula (0s|Tk|0s} is the VEV 
for the AdS bulk without boundaries and the part (T k )^ is induced by a single brane located 
at j/ = j when the second brane is absent. For j = a the single brane induced part is given by 
formula 




k\(a) _ 


k° +1 z D 6 k 


2 D - 1 n^T (§) 


K PC 

IT l 


dwu 


D-l 


Ku '(uz a ) 


Z > Z a , 


(9) 


where I u (z) and K v (z) are the Bessel modified functions. The corresponding formula for j = b 
and z < Zb is obtained from © by the replacements I K. The last term on the right of 
formula © is induced by the second brane. In the corresponding formula we use the notations 


F^[g(v)} = 


-2C 


v g (v) + [D + 


4C 


4C-1 


+ I v 2 + v 2 + 


2v 2 


D( 4C - 1) 


9 2 {v) 


vg{v)g\v)+ 


i = 0,1,..., D-l, 


D 


F^ D) [9(v)} = -—g' 2 (v) + — (4:(-l)vg(v)g'(v) + 

+ I [ w 2 + ^ + 2 (D(D + 1) - D 2 / 2] g 2 (v), 


for a given function g(v) 


G„\u, v) 


I v {v)K^{u)-K v (y)I^(u), 
Ri b \v)/Kl a \u ) 

Ki a) (u)li b) (■ v ) - lU b) ( v)li a) (u) ’ 


( 10 ) 


( 11 ) 


( 12 ) 

(13) 


and the corresponding expression for v) is obtained from formula m by the replacement 

Ku b \v)/Ki a \u) — ► Il a \u)/Ii b) (v). As we could expect from the problem symmetry the brane 
induced parts of the vacuum energy-momentum tensor corresponding to the components on 
the hyperplane parallel to the branes are proportional to the corresponding metric tensor. On 
the boundary the vacuum energy-momentum tensor for a single brane diverges, except the case 
of a conformally coupled massless scalar. The leading terms of the corresponding asymptotic 
expansion near the boundary are the same as for a plate in the Minkowski bulk. For large proper 
distances from the brane compared with the AdS curvature radius, ko\y — a| 1, the single 
brane induced energy-momentum tensor vanishes as exp[2 vkr)(a — y)] in the region y > a and 
as exp [fee (2z^ + D)(y — a)] in the region y < a. The same behavior takes place for a fixed y — a 
and large values of the parameter ko- In the large mass limit, m kp, the boundary parts are 
exponentially suppressed. For a conformally coupled massless scalar held one has v = 1/2 and 
the corresponding cylindrical functions are expressed through the elementary functions. In this 
case the corresponding VEVs can also be obtained from the results for two parallel plates in the 
Minkowski bulk by using the standard conformal transformation technique |7j. 

For the case of two-branes geometry the vacuum force acting per unit surface of the brane at 
y = j is determined by the ^ - component of the vacuum energy-momentum tensor at this point. 
The corresponding effective pressures can be presented as a sum of two terms: pj = pj \ + Pj(int )) 
j = a, b, where pj\ is the pressure for a single brane at y = j when the second brane is absent. 
This term is divergent due to the surface divergences in the vacuum expectation values and 
needs an additional renormalization. The term 

fcD+l p oo 

Pj (int) = d / n \ / dxX D ~ 1 Q. jv {xZ a /z j ,XZ b /z j ) 

2 D 7 r“T (y) Jo 

x {(x 2 - v 2 + 2m 2 /k 2 D ) B 2 - D(4£ - 1 )AjBj - A 2 ] , (14) 


3 










is the pressure induced by the presence of the second brane, and can be termed as an interaction 
force. For Dirichlet scalar these forces are always attractive. In the case of Robin boundary 
conditions the interaction forces can be either repulsive or attractive. Moreover, there is a region 
in the space of Robin parameters in which the interaction forces are repulsive for small distances 
and are attractive for large distances. This provides a possibility to stabilize interbrane distance 
by using the vacuum forces. For large distances between the branes, the vacuum interaction 
forces per unit surface are exponentially suppressed by the factor exp[2 vkr>(a — b)\ for the brane 
at y = a and by the factor exp[&t>(2zz + D)(a — b )] for the brane at y = b. In the Randall- 
Sundrum braneworld with arbitrary mass terms Cj, j = a, b , on the branes the Robin coefficients 
are expressed through these mass terms and the curvature coupling parameter as 

~ ~ 1 

Aj/Bj = --(rijCj + 4 D(k D ), n a = 1, n b = -1. (15) 

For the twisted scalar Dirichlet boundary conditions are obtained on both branes. The corre¬ 
sponding VEV of the energy-momentum tensor is defined by formulae ©, © by an additional 
factor 1/2. 


3 Surface densities 


In the previous section we have considered the bulk energy-momentum tensor. For a scalar field 
on manifolds with boundaries in addition to this part, the energy-momentum tensor contains 
a contribution located on the boundary. The surface part of the energy-momentum tensor is 
essential in considerations of the relation between local and global characteristics in the Casimir 
effect. For an arbitrary smooth boundary dM s with the inward-pointing unit normal vector n l , 
the surface part of the energy-momentum tensor is given by the formula |H. T)'/ urf ' 1 = 5(x ; dM s )rik 
with 

Tik = (<p 2 K ik - (2C - 1/2 )h ik (pn l Vnp, (16) 

where Kik is the extrinsic curvature tensor of the boundary and hik is the corresponding induced 
metric. In the problem under consideration the VEV of the surface energy-momentum tensor 
on the brane y = j is presented in the form 


(0|r«)|0> 


dtiv'Tlj 


(k D - (2C - 1/2 ) Aj/Bj 


(0|^ 2 |0), = ,,, 


(17) 


and (Olr^lO) = 0. From the point of view of physics on the brane, Eq. (1171) corresponds to the 

("suj-f) (surf) 

gravitational source of the cosmological constant type with the equation of state e) ; = —p K - . 
By making use of expressions urn it can be seen that the surface energy in the Randall-Sundrum 
braneworld vanishes for minimally and conformally coupled scalar fields with zero brane mass 
terms. The VEV (El diverges and needs some regularization. We use the method which is an 
analog of the generalized zeta function approach. The VEV of the held square is presented in 
terms of the generalized zeta function 


oo D+s (0 ! 

>- / \ 9v \^v,n^lt 'U'i/ 7 n,Zj 

0(s) = 2^ - 


i.=l Q u 9v ( uz ai uz b)\u=u l/t71 


(18) 


with j, l = a,b, l / j, gl l \u,v) = J v {v)Yf) l \u) - Y v (v)J^\u), as 


( 01 V 7 10)2=2^ — k D 


D-i^D^jPo d 


D — 1 D-l + s 


Cj( S )|s=—1, 


(19) 
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where | s= _i is understood in the sense of the analytic continuation, and B(x,y) is the Euler 
beta function. Using the Cauchy’s theorem on residues, the following integral representation for 
the zeta functions is obtained 


1 7 r f°° 

(j(s) = -sin — (D + 1 + s) / duu D+s [rijUj u {uzj ) + Bjflj u (uz a , uzb )], 

n 2 Jo 


( 20 ) 


with U au (x) = I\ v (x)/R[j l \x) and Ub u (x) = I u (x)/ffl {x). The contribution of the second term 
in the square brackets is finite at s = — 1 and vanishes in the limits z a —> 0 and Zb —► oo. 
The first term in the square brackets of this expression corresponds to the contribution of a 
single brane at z = Zj when the second brane is absent. The regularization is needed for 
this term only and the details for the corresponding analytic continuation can be found in 
[3]. For the analytic continuation of the single brane zeta functions we subtract and add to 
the integrands leading terms of the corresponding asymptotic expansions, and present them 
as sums of two parts. The first one is convergent at the physical point and can be evaluated 
numerically. In the second, asymptotic part the pole contributions are given explicitly. As a 
consequence, the single brane surface Casimir energies for separate left and right sides of the 
brane contain pole and finite contributions. The remained pole term is a characteristic feature 
for the zeta function regularization method and has been found in the calculations of the total 
Casimir energy for many cases of boundary geometries. As in the case of the total vacuum 
energy, the renormalization of these terms can be performed by using the brane tensions. For 
an infinitely thin brane taking the left and right regions together, in odd spatial dimensions 
the pole parts cancel and the surface Casimir energy is finite. In this case the total surface 
energy per unit physical volume on the brane (surface tension) does not depend on the brane 
position. The results of the corresponding numerical evaluation show that, in dependence of 
ratio of coefficients in the boundary condition, the surface energy for a single brane can be 
either negative or positive. In two branes geometry the surface energy density on the brane at 
z = Zj is presented as the sum 


(surf) 

e i 


(surf) * (surf) 


( 21 ) 


(surf) • 

where £ ' jX is the surface energy density induced on the corresponding surface of a single brane 

a t V = 3 when the second brane is absent, and the term 


Ae 


(surf) _ 


20 Bj - (4C - 1) Aj 


n 


U) ( 


(k D z j ) D B j /3 D+ 1 


r 

Jo 


duu D 1 fl 
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(uz a ,uz b ) (22) 


is the energy density induced by the presence of the second brane. Note that this part is located 
on the surface y = a + 0 for the brane at y = a and on the surface y = b — 0 for the brane 
at y = b. On the surfaces y = a — 0 and y = b + 0 the surface densities are the same as for 
single branes. The expression on the right of Eq. El is finite for all values z a < Zb and is a 
function on the interbrane distance only. In the limit of large AdS radius, kp —> 0, the result for 
two parallel plates on the Minkowski bulk is recovered. For large distances between the branes, 
the surface densities induced by the second brane are exponentially suppressed by the factor 
exp[fcc)(2iz + .D)(a — b)] for the brane at y = a and by the factor exp[2i'fcc>(a —&)] for the brane at 
y = b. The exponential suppression also takes place in the large mass limit. From the viewpoint 
of an observer living on the brane at y = j, the D-dimensional Newton’s constant is related to 
the higher-dimensional fundamental Newton’s constant by formula 


r< (-D — 2)kpGp + i (D-2)k n (b-j) 

Dj “ e (D-2)k D (b-a) _ i 


(23) 


In the orbifolded version of the model an additional factor 2 appears in the denominator of the 
expression on the right. For large interbrane separations this constant is exponentially small on 
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the brane y = b and the gravitational interactions on this brane are exponentially suppressed. 
The corresponding effective cosmological constant generated by the second brane is determined 
by the relation 


A Dj = 8*G Dj Aef m ' f) 


Kf 


(24) 


where Mjyj is the D-dimensional effective Planck mass scale for an observer on the brane. For 
large interbrane distances the quantity (PI is suppressed to compared with the corresponding 
Planck scale quantity in the brane universe by the factor exp[fc£>(2^ + D)(a — b)\, assuming that 
the AdS inverse radius and the fundamental Planck mass are of the same order. In the original 
Randall-Sundrum model with D = 4, for a scalar field with the mass |m 2 | < kp, and interbrane 
distances solving the hierarchy problem, the value of the cosmological constant on the visible 
brane by order of magnitude is in agreement with the value suggested by current cosmological 
observations without an additional fine tuning of the parameters. 
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